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3^ ! Abstract 

This paper serves as a report on the large amount of analysis done in conjunction with the 

(N 

^ ' polarized proton program at RHIC. This comprises elastic scattering data of protons on protons in 

cn 

0^ ■ colliding beam or fixed target mode and proton beams on carbon targets. In addition to providing 

in 

'^ I a model for the energy dependence of the analyzing power of elastic scattering needed for proton 
polarimetry, it also provides some significant information about the spin dependence of dominant 

^D ■ Regge poles. Most notably, the data indicates that the Pomeron has a significant spin-flip coupling. 

- ^^ , This allows the exploration of the double spin flip asymmetry A^r^ for which some data over a 

X 

?H ' wide energy range is now available, along with a concrete realization of a proposed Odderon search. 
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1. INTRODUCTION 

Recent experiments carried out at RHIC using polarized proton beams give us the oppor- 
tunity to address a question of long-standing interest; namely, what is the spin-dependence 
of the Regge couplings to the proton and, in particular, does the Pomeron coupling depend 
on the helicity of the proton? The work described in this paper results from a program of sev- 
eral years duration; background information and earlier results can be found in [l|, y, l3|, y, |5| . 
The program originated in the attempt to calculate the energy dependence of the analyzing 
power A]\}{t) — the single transverse spin asymmetry — for the elastic scattering of polarized 
protons off protons or nuclei, especially carbon. This analyzing power can then be used for 
proton polarimetry in a high energy colliding beam facility for energies where no calibration 
standard exists or where rapid measurements are required. The goal of the present paper is 
to describe the physics that has been learned from this work. 

Starting from the well-justified assumption that high energy elastic proton-proton scat- 
tering can be described accurately with five Regge poles, the Pomeron and the dominant 
C = ±1 poles for / = 0, 1 [6], we will use Ajs[{t) measured at two energies in pC scatter- 
ing and at one energy in pp scattering to determine the spin dependence of the five Regge 
residues in this description. As we will see, the analysis with this limited number of pa- 
rameters requires a knowledge of the beam polarization P at only one energy; at the other 
energy in pC scattering, only the shape of the analyzing power in t in the Coulomb Nuclear 
Interference region (CNI) is required. The method makes use of the spin-dependent asym- 
metry Ai\f{t) in the elastic scattering of two nucleons induced by the interference between 
the 1-photon exchange amplitude and the strong, hadronic scattering amplitude. The sin- 



gularity of the so-cal 



ed Coulomb amplitude leads to a characteristic enhancement of A^lt) 



for very small —t [l|, [Zl, |8|, |9(] . This enhancement is important in yielding a practical signal. 
In order to separate the 1 = and 7 = 1 exchanges for the Regge analysis it is necessary 
to have two targets; fortunately we have data for proton and carbon targets. Key to using 
pp and pC measurements together to extract physics (as opposed to simple polarimetry) is 



the relation, proved to be valid under a wide range of assumptions [10(], that the ratio of 
single-flip to nonfiip amplitudes, commonly called r (see Eq. 2.14 below), is the same for 
both processes; more precisely Tpc = Tq ~ {rpp + r„p)/2. 

Because these calculations were going on simultaneously with the experiments they are 



designed to describe, various ways of using the underlying Regge theory in order to extract 
the needed information from the hmited data were tried [2|, l3|, U, |5|. Some of these will be 
described here because they shed some light on the validity of the as-yet-untested theory we 
are using; furthermore they may prove useful in other circumstances in the future. 

This paper is organized in the following way. Section 2 will be a brief review of the general 
CNI theory used here, including quasi-elastic nuclear scattering. In Section 3.1 we will apply 
it to pp data obtained at RHIC with proton beams at 24 GeV/c and 100 GeV/c on fixed 



proton targets 
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12|. Errors will be estimated for the spin-flip factor t{s). In Section 3.2 



we will do the same thing for pC RHIC (and AGS) data, which came chronologically mostly 



before the pp data was obtained [13|, llj|. Because the target is a nucleus, the analysis is 
somewhat more complex than the pp case and we will describe that here. Errors will be 
given here, too. In Section 4 we introduce the energy dependence using the Regge model 
of |6| and determine the energy dependence of the pp flip factor r(s) and the pC flip factor 
To{s) using several different methods. All of the Regge poles, in particular the Pomeron, 
are found to have significant spin-flip coupling. The pC fit gives very good determination 
of the 1 = Regge couplings, including the Pomeron. The pp couplings are less certain. In 
Section 5 we will assume that the fit coefficients in [6| are given by factorizable Regge poles 
and determine the spin-flip couplings of the p, u, f, and 02 . In Section 6 we will turn to the 
new and limited data on the double-spin asymmetry Aj^f^ [iSl, ll6| and see what our model 
has to say about that. It will be necessary to take account of Regge cuts since factorization 
of pole couplings forces their contribution to Ann to vanish at t = 0. This leads naturally 
to an exploration in Section 6.1 of the idea presented by Leader and the author ^] for 
searching for the Odderon. Finally Section 7 sums up. 

2. GENERAL CNI DISCUSSION 

The method used here has a long history dating back at least to [7]. In that paper, 
Schwinger proposed using the same effect to produce a beam of polarized neutrons. {Almost 
the same. The neutron singularity at t = is milder because it is neutral, but the magnetic 
moment provides the needed enhancement.) The details for pp and pC are slightly different 
and we summarize the needed formulas below. For an overview of polarization phenomena 
see Ui 



2-1. Proton-proton CNI 



a 



A detailed discussion can be found in [1]; a summary follows. Five independent helicity 
amplitudes are required to describe proton-proton elastic scattering [9|, [l9| : 
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(2.1) 



Here we use the normalization of 91]. Since we are interested only in very high energy ^/s, 
such as will be available at RHIC, and very small momentum transfer |t| < 0.05 GeV^, we 
will generally neglect m with respect to s and neglect t with respect to m to simplify the 
presentation of the formulas which follow. The total and differential cross sections are given 

by 

An 
atot = —ImiMs,t) + Ms,t))\t=o (2.2) 

s 

and 
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(2.3) 



Using only initial state polarization, with one or both beams polarized, one can measure 
seven spin dependent asymmetries. We follow the notation of |9| . There are slight variations 
in the definitions used in the literature, having to do with the orientation of axes. 

^7v— = ^Irn{(j)l{(j)i + (f)2 + (j)3-(t>i)}, 
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Ann^ = ^{2|05r + i?e(0t02-03</'4)}, 
Ass^ = ^i?e{0i0; + 0304}, 
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It will be convenient to introduce some shorthand: 



^n 



(2.4) 
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(2.5) 



which enter into the two cross section differences corresponding to longitudinal and trans 
verse polarization: 

Im0_(s,O) lAcx^is) 



Im0+(s,O) 2 crtot(s) 

Im02('5, 0) Acr^(s) 



Aa^ = a^-a^, (2.6) 



Im0+(s, 0) o-tot(s) 

At these small values of t, the interference of the strong amplitudes with the single photon 

exchange amplitudes will be important; this interference is central to this paper. To lowest 

order in a, the fine structure constant, one replaces 

0, = 0^^^ + <PreM^S)GE{q'f (2.8) 

with hadronic and electromagnetic elements. The Coulomb or "Bethe" phase S is approxi- 
mately independent of helicity [ol, |20| 

^ = "^ %^(E + 8/A^) -"^ ^2.9) 

where B is the logarithmic derivative of the differential cross section at t = and we use 
the fit 

S(s) = ll + 0.51n(s/102) (2.10) 

in GeV~^ through the RHIC region. Also g^ = —t, Euler's constant 7 = 0.5772 . . . and 
A^ = 0.71 GeV^ reproduces the small momentum transfer dependence of the proton form 
factors assumed to satisfy 

GE{q') = GM{q')/f^p = (1 + qVAY'- (2-11) 

For pp scattering at high s and small t, the electromagnetic amplitudes are approximately 

0r = 0r = y, 

^em _ -em _ ^^K^ 



2my/—t 
where Hp = k+1 is the proton's magnetic moment, and m its mass. For the full expressions 



Is 



see, e.g., [9|. For our purposes it will suffice to write 



0+(.,t) + 0_(s,t) = ^a,ot(5)(^ + p(3))e^(^)*/^ (2.13) 



where over the tiny CNI t-range we will neglect any variation of p. Likewise with 



b,{s,t) = T{s)^^Ms,t) (2.14) 

m 

we will assume that r, which is complex, varies with s but not with t over the CNI range. 

(Sometimes in the literature one uses r^{s) = (l)^{s, t)/{^J —t/im? Im0+) = r(s)(p(s) +i). We 

prefer to use r.) Then for simplicity for pp scattering if we neglect 5 and the amplitudes 

4>2 ,4'A, which are not enhanced by the Coulomb singularity and are probably small (see 

Section 6 below), the equation for A^ in Eq. (2.4) gives 



Snm da^_j,^ANis,t) _.,^ ^_,, ..,_ ...tc ,t™ .^,, , .^2^ 



e-i^tl^lipiL = (k/2 - Rer(s) - p(s)Imr(s))- + Imr(s)(l + p(s)^). (2.15) 



O-tot(s)^ dt ^/^ t 

tc = —Sira/atot- This formula will be used with the recent data for AN{s,t) to determine 
the (complex) spin-flip factor for pp scattering t{s). 

2-2. Proton-carbon CNI 

The basic physics here is the same as for pp, but it is simpler to describe because there 
are only two ampltudes For pC scattering the equation has the same form with all the 
quantities, including tc modified to refer to the carbon charge and wave function. Because 
the carbon is rather large, there will be significant variation with t for ppc and for the ratio 
of the electromagnetic to hadronic form factors. Also the Bethe phase is significant and 
needs to be kept. This was all worked out in 1 101] using the harmonic oscillator carbon wave 



function used by Glauber and Matthiae 21 
single particle densities 



22l. The needed functions are: 
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(2.16) 



the corresponding electromagnetic form factor 

F'=™(t) =Aac [ dbMV^b)be'''^''\l + \ ac{b^ + — )) 
J 3 ac 



(2.17) 



the hadronic amplitude 

FHit) = ImFo^(g)/ImFo^(0) (2.18) 



where for carbon 23| 



F^{q) = % \ d'he"'-" { 1 



fihJ^-^ 
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(2.19) 



Note that p in this expression denotes the / = nucleon-nucleon real-to-imaginary ratio, 
very close to the pp value, ac = 0.0152 GeV^. The analog of Eq. 12.151 is 



with 



where 



with 



16 TT d (TpA 



[cr 



i1y dt 



<{t) = ^ F\{t) iFrit) \ ^{l- S,A P,a) 



- 2 Re To - {6pA + Ppa) )] - 2 F^{t) (im ro(l + p^^) ) | . 
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tc = -8TrZa/a^„^. 



PpAis^t) = ReF^iq)/ImF^iq). 



(2.20) 



(2.21) 



(2.22) 



(2.23) 



1. Inelastic corrections 



The above calculation was used in the paper with Kopeliovich lOj and all of my reports 
prior to 2004 [2|, y, |j]. Just after the last of these a preliminary measurement of the pC 
differential cross section was reported. Although it is still preliminary, it is clear from the 
data that the differential cross section calculated in 10| does not agree with it. This becomes 
very clear above the CNI region where at diffraction dip predicted at t = —0.09 becomes 
in the data a break in the exponential fall which has a slope of about 60 GeV^^ between 



-.02 and -.06 241 . Exactly this effect was reported long ago in the CERN proton-nucleus 
measurements |25j where it was explained by additional incoherent proton-nucleon scattering 
which comes in because of the resolution. Subsequently, Glauber and Matthiae |22| extended 
the Glauber method used by 10| to take into account quasi-elastic scatterings. It is based on 
completeness and only requires that the experimental resolution be wide enough to include 
nearly all excited nuclear states. The mass resolution presented at Blois 2005 by Bravar is 
about 1 GeV 26|. Since the binding energy of carbon is about 92 MeV all of the excited 
states that do not lead to production or breakup must be included and the completeness 
assumption is good. We will expand in powers of the scattering and take just the first term 
to estimate this correction. The key parameter is Ni{s) which determines the number of 
quasi-elastic scattering. According to [22] it is given by 



NAs) 



2tx 



(2.24) 



and the nuclear thickness function T(6) is given by the integral over the nuclear density 



T(6) = 4( 



OCn 3 



TT 






(2.25) 



A'^i(s) varies slowly with s and is approximately three quasi-elastic scatterings throughout 
the energy range. The measured pC cross section is then the sum of elastic cross section 
Eq. 12.211 plus A'^i(s) times Eq. l2.3[ The comparison between the elastic cross section and this 
sum is shown in Fig. 1. The difference is not large, but it is significant at the largest —t. 
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FIG. 1: In-^ at p^ = 24GeV/c. The lower curve is from [lOJ], the upper curve is based on [22l | and 



includes quasi-elastic corrections 



For completeness Fig. 2 shows a larger range including the diffraction dip. An exponential 
fit to this calculation between -.02 and -.06 gives a slope of 58.2 GeV~^, close to the CERN 
slope |25l |. 
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FIG. 2: In -^ at p^ = lOOGeV/c over a larger range of t. The lower curve is from [10|, the upper 
curve is based on ^] and includes quasi-elastic corrections 



SPIN FLIP FACTORS 



3-1. pp 



Ai\[{s, t) has been measured at RHIC for proton beam energies of 24 GeV/c and 100 GeV/c 



on a polarized gas jet target 
the results for r^ reported 

r(s) =r5{s)/{i + p{s)): 



271]. The data has been fit to the formula given in Section 2 and 



u, 



12] • We have converted these to our preferred parameter 



r(s]100]) 
r(s[24]) 



-0.0148 + 0.002i, 
-0.100 + 0.0306Z 



(3.1) 
(3.2) 



s[p] = 2m? + 2m ^Jp^ + m?. (Because it occurs so often, we will often write r(p) instead 
of r(s[p]) when there is no chance of confusion.) The quoted chi-squares for these deter- 
minations are: at 100 GeV =11.1 for 14 points and at 24 GeV = 2.87 for 9 points. The 
error matrices were also given and used in our calculation of the la error ellipses, shown 
in Fig. 3 and Fig.4. The best fit at 100 GeV/c has r consistent with 0, as reported by the 
experimentalists; the error ellipse for the much shorter run at 24GeV/c is about 3 times 
larger than the 100 GeV/c case but the value of r is significantly non-zero. This was the 
last data to become available and had a decisive influence on our results. The errors will 
surely be improved in the future. 
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Re T(IOO) 



FIG. 3: lo" error ellipse for the proton r determined with pi = 100 GeV/c beam on a polarized 
gas-jet target. 

Im X(24) 
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FIG. 4: Icr error ellipse for the proton r determined with p^ = 24GeV/c beam on a polarized 
gas-jet target. 

3-2. pC 

This process has a longer and more complex history than that of pp. The first results 
come from an AGS experiment E950 [13| at 21.7GeV/c. There were also measurements at 
RHIC at 24 GeV/c and 100 GeV/c. The beam in the E950 experiment had its polarization 
determined by a comple tely independent experiment E925 based on a known analyzing 
power for pp scattering 28|]. The value they used was 0.407 ± 0.036(stat) ± 0.049(s?/st). 
For a few years the RHIC measurements above injection energy used only estimates until 
2004 when the polarized jet target became available at 100 GeV/c. At that point a really 
extensive set of data with very small errors was obtained. In general, the experimentalists 
analyzed their data using Eq J2.20l and reported r values. (In fact, for historical reasons, they 
usually report values for r^pc = To(PpC + 0- Here we will always convert those numbers to 



9 



r.) Because we now know that there are corrections due to quasi-elastic scattering, we have 
reanalyzed their data using the corrected form for the asymmetry: i.e. for each polarization 
state of the proton beam we calculate the sum of the corresponding elastic pC differential 
cross section plus the A^i times the pp differential cross section. The difference of these two 
divided by the sum is the corrected analyzing power as a function of tq. The value of Tq 
determined in this way is a little different from the uncorrected fit from E950 [l3|, but the 
chi-square is very much better, 0.55/d.o.f. compared to 1.16/d.o.f . The same comparison 
for the 100 GeV/c data gives 1.32/d.o.f. for the corrected fit compared to 1.51/d.o.f for the 
uncorrected fit. We will use the corrected fits for each of the data sets we use. The results 
for E950 and for RHIC at 100 GeV/c are 

ro(21.7) = -0.222 - 0.0584i (3.3) 

ro(lOO) = -.011 - .0498^ (3.4) 

The subscript indicates that this is the r value for / = exchanges. (At 100 GeV/c we 
have used only the lowest 10 values of — t, omitting the last 4 points they report. The chi- 
square is considerably better that way and the corrections to the low t approximations used 
in deriving Eq. 12.201 are smaller.) The error ellipses are shown in In Fig. 5 and Fig. 6. From 

ImTo(21.7) 




-0.4 -0.3 -0.2 -0.1 



ReXo(21.7) 



FIG. 5: pC error ellipse for E950 

these figures we easily see that the spin-flip factor is non-zero at both energies. Indeed, at 
21.7 GeV/c the flip factor is bigger than 20%, and there is significant energy dependence. 
The fits to the two data sets are shown in Fig. 7 and Fig. 8. (The data for E950 are from 
[13|. The data at 100 GeV are from [M])- Notice how precise the latter set is. As a result 
it will play a central role in our parameter determination. The chi-square for 100 GeV fit is 
10.6 for 10 data points, significantly better than the best fit without quasi-elastic corrections 
which is 12.1 for 10 data points. 
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FIG. 6: pC error ellipse for 100 GeV/c 
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FIG. 7: pC analyzing power for E950 
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4. ENERGY DEPENDENCE OF THE SPIN-FLIP FACTORS t(s) 



It would be very useful to know the energy dependence of r(s) so that the analyzing 
power could be used for polarimetry at energies where it has not been measured. At a 
deeper theoretical level, it should giv e us some insight into the old problem of the spin- 



dependence of Regge couphngs [29|, l30|, l311- We propose to take here a very simple approach 



which is appropriate for the e nerg y and t range of the physics under consideration. For an 
overview of Regge theory see 32|]. From among all the various Regge fits, we have chosen 
to use the one given by the Particle Data Group due to Cudell et al. \&\. This is based on 
a Regge picture with two simple Regge poles with C = ±1 with 0<J<latt = plus a 



An 
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FIG. 8: pC analyzing power for pL = lOOGeV/c with data from 
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more complex structure at J = 1 representing the Pomeron. (In the earher reports 2|, y, |4| 
we used an alternative fit from [6] in which the Pomeron was parametrized as a simple pole. 
This contributes to the difference between the spin-flip factors given here and those in the 
earlier reports.) The formula for the forward pp amplitude 

— 0+(s, 0) = gpis) + g+is) + g^s) (4.1) 

= 9o{s) 

with 

gp{s) = 0.789 TT In (s/29.1) (4.2) 

+i(91.26 + 0.789 ln2(s/29.1)), 
g+{s) = -Y,s-^{cotn/2{l-r])-z), 
g4s) = -Y2S-"'itann/2il-r]')+i), 

with Yi = 109.51 GeV"^ Fa = 85.86 GeV^l The parameters are rj = 0.458, r/' = 0.545 and 
the coefficients are all in units of GeV~^. g^ and g^ arise as sums of contributions of C = +1 
(/ and a2) and C = — 1 (u; and p) of / = and 1 = 1, respectively, assumed in each case to 
be degenerate in J. 

To introduce spin dependence we write 



— 05(s,t) = ^-^{rp(7p(s) + r+(7+(s) +r_(7_(s)} (4.3) 

O lib 

For simple Regge poles tr is real and we will assume the same is true for the Pomeron. Then 

^(A - 'Tp gpjs) + T+g+{s) + T^g-{s) 
^' gp{s) + g^s) + 9-{s) ^'^ 

and the energy dependence of the real and imaginary parts of r are given by three real 

parameters. 

The pC scattering is pure / = and so requires a different set of functions: presumably, 

the Pomeron amplitudes are the same but the C = +1 amplitude is pure / and the C = — 1 

amplitude is pure u so g+ and g_ are replaced in the pp amplitude to get the pC amplitude 

by gf and g^^, respectively. From [6|] 

gf{s) = -Yf s-^{cot 7r/2(l - r^) - z) (4.5) 

gM) = -Y^s~^'itann/2il-ri')+i). 

12 



with Yf = 106.36 GeV"^ and Y^ = 81.49 GeV"^ not so different from Eq. EM because the 
spin-independent parts of pp and pn scattering are nearly the same. 
In obvious notation 

4-1. Determination of the 1 = spin-flip factors 

The measurement of the analyzing power at any energy s determines two parameters, 
Rer(s) and Imr(s) and so knowledge of An at two energies is more than enough to fix 
the individual spin-flip factors. Early in this work we encountered a situation where the 
asymmetry was measured but the polarization was not known. By separating the singular 
and non-singular terms in Eq. 12.201 one can determine the shape of the asymmetry So{s), 
in t. This is independent of the polarization P and provides another relation between Rer 
and Imr at that energy: 

Ms) « p,"'"f\,y (4.7) 

P{k/2 — Rero(s)) 

Therefore knowledge of the analyzing power at one energy and the shape at another is 
sufficient to determine the flip factors and, thereby, the polarization/analyzing power at the 
second energy. 

1. E950 T and the shape of the first 100 GeV RHIC data: a test 

Historically, the flrst application of this method was to the flrst reported data (2002) for 



polarized protons on carbon at RHIC J33j. In addition to data with a proton beam whose 
polarization was measured using the E950 calibration, the asymmetry was also measured at 
100 GeV where the polarization was unknown. We flt data in Table 1 taken from [33], to 
the formula in Eq.f l2.20l) with the right hand side multiplied by the unknown P(IOO). The 
result of the regression is 

P(100)(l--Rer(100)) = 0.263 (4.8) 

K 

P(100)-Imr(100) = -0.0137. 

Combining these together we get for the shape of the distribution 

01 37 
^o(lOO) = -^^ = -0.052 (4.9) 
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-t 


e{t) 


e{t) 


0.0117 


0.0036 


0.00055 


0.0138 


0.0029 


0.00047 


0.015 


0.0034 


0.00060 


0.0184 


0.0018 


0.00069 


0.0194 


0.0014 


0.00058 


0.0217 


0.0025 


0.00087 


0.0249 


-0.0004 


0.00072 


0.0251 


0.0009 


0.00085 


0.0306 


-0.0010 


0.00072 


0.0360 


0.0010 


0.00091 


0.0416 


0.0013 


0.00116 


0.0473 


-0.003 


0.00146 









TABLE I: -t, raw asymmetry e{t), and errors e{t) for RHIC 100 GeV/c 

By expesssing the value of r(21.7) as determined in E950 in terms of the Regge spin-flip 
couplings, and in the same way express 50(100) in terms of the Regge spin-flip couplings 
via ro(lOO), we have three equations to solve with the result Tptcst = 0.028 ± 0.14, rjtest = 
—0.967 ±0.35, Tf^tcst = 0.509 ±0.23. There are significant errors in these determinations and 
we will move on to a more accurate determination in the next section. 

We could now calculate rotest(s) at any higher energy, with corresponding accuracy. (The 
model as it stands, is not really suitable for going to lower energy because lower lying 



Regge poles wi 
this limitation 



1 rapid 



30, 



M, 



z. 



beome important. Thanks to Boris Kopeliovich for emphasizing 



35|.) 



To test and illustrate the model, we will use these results first for rotest(lOO) and thereby 
determine the polarization at 100 GeV/c, P(IOO). Using these values of the Regge coupling 
we get ro(lOO) = -0.097±0.15-(0.054±0.02)i. With errors this large this method will give 
only a crude indication of P. We find P(IOO) = 0.23 ± 0.08. This is a little smaller, about 



15%, than the value 0.27 used by [331] based on assuming An to be energy independent with 



the E950 value. In Fig. 9 we show the raw asymmetry measured at 100 GeV/c plotted with 
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the prediction using Eq. 2.19 and P = 0.23. The agreement is reasonable, with x^/dof about 
1.5. 




■ t GeV2 



FIG. 9: Raw asymmetry e at 100 GeV/c l3S] with curve predicted by the model ro(lOO) and polar- 
ization P, predicted to he 0.23 



2. Shape of E 950 and precision measurement at 100 GeV 



The errors on P and r determined in the previous section are necessarily large because 
the errors on both data sets are large. We should be able to do much better by using 100 
GeV T values determined in Section 3.2 which have very small errors and use just the shape 
from E950, which has a much smaller error than r0(21.7) itself. The measured value of the 
shape is S'(21.7) = —0.0523, using the quasielastic corrected Ajy, and the error on it can be 
determined from the propagation of errors from the regression of the 21.7 GeV data. It is 
very small: (5Shape(21.7) = 0.003. 

We can use that along with the error matrix for r(lOO) to determine the error matrices 
and ellipses for the three Regge spin-flip couplings. Fig. [TOl The errors are seen to range 
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FIG. 10: Regge coupling error ellipses 

from about 6% to 13%. 

We have tried alternative determinations of Tp using the 100 GeV/c data and the 21.7 
GeV/c data in various ways with this model, but this one has the smallest errors, and it is 
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certainly non-zero. We will use the determination of residues here to calculate tq{s) in the 
following. The values are 



Tp = 0.10 ± 0.01, Tf = -0.79 ± 0.05, r^ = 0.52 ± 0.06. 



(4.10) 



Note that here Tp is clearly non-zero, while in the "test" case it is consistent with zero (as 
well as with this non-zero value.) 

These give the energy dependence of Tq as shown in Fig. 11. 
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FIG. 11: The energy dependence o/Rero andlrriTQ through RHIC fixed target range with la errors. 

We can use the Regge coupling error matrices to determine the errors on the prediction 
for r0(21.7). Fig. 12. The central value is consistent with the experimental measurement 
which has, for the real part, a much larger error than that of the prediction. 

There also exists a set of pC asymmetry e data taken at 24 GeV/c at RHIC, but the 
final results (with errors) have not been released and so have not been used in the present 
analysis. 

4-2. Determination of the proton spin-flip factors 

We will describe here two distinct ways of determining the proton spin-flip factors. If the 
model and the data were perfect, they should give the same results. 



1. 100 GeV/c pp measurement in conjunction with Tp from pC 

The first pp data that was available to us was the 100 GeV/c data reported in Section 2. 
That provides us with two new parameters to determine the two C = ±1 spin-flip residues 
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FIG. 12: Error ellipse for predicted r(21.7) shown within the measured ellipse as determined by 
E950. See Fig.5 

which we label r_|_ and r_. Assuming the model is correct, the same value of Tp applies for 
the pp case and so we have enough information to determine r(s) for pp. 



Tp = 0.10, r+ = -0.51, r_ = 1.16 



(4.11) 



with the error ellipse for the two new couplings shown in Fig. 13. These values are determined 




T+ 



FIG. 13: Error ellipse for pp Regge couplings with C = ±1. 

by the spin-flip couplings of / -|- 02 and oi uj + p respectively in a way to be made precise 
shortly. Because these are very different from tj and r^ we see that the 02 and p couplings 
must be very large. This has been noticed from earlier data at lower energy and from vr^A^ 
data [30,]. We also see that these new couphngs are not very well determined in spite of the 
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very nice pp data. This will surely improve in the future. On the other hand, the implied 
statistical error AA^ on the predicted A^ is not so bad. See Fig. 14. 




0.005 0.01 0.015 0.02 0.025 0.03 - t GeV^ 



FIG. 14: The la error on A^- implied by the errors on r as determined in this section. 

We used these parameters to predict A^ at 24 GeV before it was measured. The com- 
parison to data is shown in Fig. 14: This shows the model is not bad: good enough to use 
for 10% polarimetry but not for 5%. 

2. Joint fit to T 's for pp at 24 and 100 GeV/c 

We use the measured values with known polarization to determine the three unknown pp 
spin-flip factors, and determine the r's that minimize the sum of the chi-square functions 
using the measured correlation errors. The results of this procedure is 



Tp = .068 ± 0.054, r+ = -0.444 ± 0.443, r_ = 0.897 ± 0.611 



(4.12) 



These values are consistent with the values using the previous method, but the errors are 
very large. We recall that the most accurate value of Tp = 0.10 ± 0.01 was obtained by the 
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FIG. 15: Data for Aj\f at pi = 24GeV/c and prediction using parameters determined from 
100 GeV/c pp measurement in conjunction with Tp from pC. Dashed lines are Icr errors. 



pC analysis. Some have imagined that t{s) = 0; this is consistent with the 100 GeV data. 



but at 24 GeV it has a chi-square of 35.5 15|. It is difficuh to do better than this with 



the small values of r with large uncertainties. Other attempts were made: for example by 
fitting the shape of the pp data at 24 GeV with r(lOO) together or by jointly fitting pp and 
pC. Somewhat different values of the Regge spin-flip factors are found but always within 
the errors which were found to be even larger. 

4-3. Predictions for higher energy 

1. 250 GeV 

It is expected that there will be a run with polarized protons at 250 GeV/c in the near 
future so we present here the prediction using the value for r given in Section 4.2.1. Although 
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FIG. 16: Prediction for A^ at pi = 250GeV/c (long dashes) compared to fit at 100 GeV/c (short 
dashes) 

r remains small at 250 GeV/c, the change from 100 GeV/c is relatively large. The predicted 
energy dependence of r over the RHIC flxed target range is shown in Fig. 17. 

2. pp2pp at s = 200^ . 

There will soon be runs with colliding beams of polarized protons at various energies. The 
flrst one has already taken place with 100 GeV/c beams in both rings. The prediction of our 
model with the parameters just determined is t{200'^) = (0.093±0.012) + z(-0.002±0.0008). 

In 2005 the results of the flrst measurement of Aj\f in the collider mode at RHIC were 
reported 36|. The spin flip parameter r^ was determined with large errors to be r^ = 
(—0.033 ±0.035) + z(— 0.43 ±0.56) with the correlation being such that a zero value for both 
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FIG. 17: Energy dependence of Rer and Inir with errors through the RHIC fixed target energy 
range. 

lies within the la elhpse. To compare with our prediction we convert our prediction for r 
to a prediction for r^ using the calculated ppp (200^) = 0.127 : 

r^prediction = 0.0138 ± 0.008 + z(0.093 ± 0.012) (4.13) 

The comparison is shown in Fig. 18 where we have calculated the error ellipse for r^ from 



361]; the prediction, shown as a black ball to indicate approximate errors, is a little bit 



outside the la ellipse. 



Re r5 (lOO^) 




Imr5 (2002) 



FIG. 18: Comparison of the value of r^ with errors from pp2pp with the model prediction with 
estimated errors marked by a black ball. 

This is a very big step in energy, and so we are happy that the measurement and our 
model prediction are compatible. However, given the small errors of the model's predicton, 
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more accurate A^ measurement is needed at high energy. This is expected to be achieved 
in the near future. 

5. REGGE COUPLINGS 

We would hke to use these results to calculate Aj^^, for example. In order to predict 
A^j^ we need to extract the factorized residues. There is not quite enough data to do this. 
The pomeron coupling is completely fixed but the Fi term in g{s) gives i3f + l3l2^ t+^i gives 
Wi^F + PLPa2), 'TfYf gives PfPf- Similar for C = -1 poles: Y2 gives Pp + Pi, t_Y2 give 
i.PlPp'^ PLP^) ^'^'^ ''"t^K; gives /3£/5^. We still need to pull out the p and a2 non-fiip couplings. 
We use {Y2 — Y^)/2 to give us /^^ and (Yi — Yf)/2 to give us /5^2 so everything is determined. 
(Since there is not nearly as much data for pn and the energy range is much smaller than 
for J9J9, there must be some hidden errors in this parametrization). 

What to use for Pomeron? It doesn't have standard form here, but it doesn't matter for 
our program because its fiip-factor Tp is determined from the fit to pC between p = 24 and 
p = 100. We will assume that, as with a factorized Regge pole, that the double flip pomeron 
will have Tp{—t/m'^) before absorption. 

In this section we will use the central values that were given (with their errors) in the 
preceding sections. 

(5.1) 



Tp = 


= 0.10, 


^.f = 


= -0.79, 


Tu^ = 


= 0.516, 


r+ = 


= -0.510 


r_ = 


= 1.163. 



We use these values plus Yi,Y2,Yf,Y^ to produce the table 

There are a few points to be rnade here: although we have not attempted to use the same 
normalization as earlier works 
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3l| we can see similar patterns here: the / and 00 non-flip 
couplings are approximately the same as are the p and a2; the / = non-flip are much 
larger than the 1=1 but the 1 = 1 flip are very much larger than the 1 = flip couplings. 
Unfortunately, the errors here for 1=1 flip are very large. There are two reasons for this 
large error: (1) the ratio for flip to non-flip for 1 = 1 has fairly large errors, about 30%, see 
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coupling 


value 


error 


/?/ 


10.312 


0.100 


Pu 


9.027 


0.035 


Pa2 


1.774 


0.582 


Pp 


2.09 


0.150 


Pi 


-8.147 


0.518 


pL 


4.676 


0.591 


pL 


15.877 


12.729 


Pi 


27.652 


13.424 



TABLE II: The non-flip and flip residues for the pure Regge pole model as normalized in the section 

Fig. 13, and (2) both 1 = 1 non-flip couplings are rather small amplifying these errors in the 
calculation using factorization of the spin-flip factors into flip times non-flip residues. These 
errors will undermine our attempts to calculate the double-spin flip amplitude AjvAf which 
is based on this factorization. 

6. DOUBLE TRANSVERSE SPIN ASYMMETRY 



The natural next step to take in applying this model is to see what it says about A^m- 
During the recent RHIC data taking some measurements of this quantity have been made 
in pp scattering at three energies: 24 GeV/c and 100 GeV/c protons on a fixed hydrogen 
gas jet target 15|] and in colliding 100 GeV beams 16|. The asymmetry is calculated using 
Eq. 12.41 Note that near t = A^vAf and Ass become nearly equal and in 16| they measure 
both. 



We will write in analogy to Eq. 14.31 and [30|, ImI] 



1 



i2(s,t) = iz{-rlsgP{s) + {py+Pii) 



sinvra/ '^ '' sinvra 



(6.1) 



-t) 



m^ 



This is the form given by the pure Regge pole model. The physics will require that this 
be modified soon. First, parity conservation of the hadronic interaction and factorization 

n 

of the poles requires that (f)2{s,t) = — 04(s, t) [1] . At the same time angular momentum 
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conservation requires that 04(s,t) vanishes with t as t — > 0. This forces 02(s,t) to vanish 
in the same way and so kinematically suppress A^n. Because the poles are suppressed, 
they cannot be the leading contribution near t = and so we must include Regge cut 
contributions. Because the cuts break factorization, the equality of 02 and —04 is broken 
and the suppression is lost. 

Of course, bringing in cuts takes the model dependence of this calculation to a new level, 
so we will adopt the simplest possible approach to estimating cut effects, the "absorptive 



cut" model of [37|| . This has been widely used in unpolarized calculations with some success, 
so we can hope that a reasonable estimate is obtained in this way, but claim nothing more 
than that. 

We will here follow the general ideas of Kane and Seidl [3^] for absorption corrections 
to Regge poles. There is assumed to be one absorption factor S{b) in 6-space (impact 
parameter) for all amplitudes. It will be small for 6 = and rise to near 1 as 6 —i> 10 GeV~^. 
It is derived from elastic rescattering plus an artfully constructed contribution from inelastic 
scattering. I will not try to model that but simply construct a simple form for 5* that looks 
similar to what they calculate. We will need to correct most importantly the double-flip 
amplitude, but we will also need to calulate the effect on non-flip and single-flip. This will 
change the relation between the Cudell parameters and t(s), and the Regge residues. We 
will continue to assume that all terms have the same e'^'-'^^*/^ behaviour. 

We will assume here that 5* is real to avoid too many parameters. So deflne 

SKib/Ra) = l~Ke'^^ (6.2) 

i^ is a real number that determines the strength of the absorption, and we will take R^ = 
2B{s). We will start from (pi and 05, transform the simple pole forms to 6-space, then 
convolute with S. This will change the relation between the Regge residues and the Cudell 
parameters as determined in 0i and 05. We treat the Pomeron slightly differently from the 
other Regge poles: since we began with a Pomeron which is manifestly not a Regge pole, 
we assume that it already has the absorption cuts taken into account in both 0i and 05. In 
addition we assume that its spin dependence factorizes so the form is exactly as in Eq. 16.11 
with the Tp as determined from the An flts. We will plug these new parameters into 02 and 
then run it through the absorption machinary. The most important result of this will be to 
replace the —t factor in 02 by a constant times 1/B{s). 
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We have used two different values oi K: K = 1 corresponding to complete absorption 



and K = 0.6, the value favored by 38|. K could in principle depend on s but the level of 
accuracy of our other parameters and of the current data does not really allow a serious 
consideration of this. The result of this simple model is that now the non-flip residues are 
increased by a factor of l/v^cTS for i^ = 1 or l/v^O/f for K = 0.6 from the values determined 
by the Cudell parameters. At the same time the spin-flip residues are decreased by a factor 
of 0.94 or 0.98 for K = 1 or K = 0.6, respectively. The most important change is the 
behaviour of 02 as ^ ^ 0; it now no longer vanishes but, in the small t region, —t is replaced 
by j^-^ + 0{t). The same calculation yields 



<l>T\s,t) = B{s)t/8<PTV,z 



(6.3) 



explicitly breaking the factorization. The model now predicts a significant value for A 
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FIG. 19: Model prediction of Ann for 24 GeV/c and 100 GeV/c protons on gas jet target for 



absorption parameter K = 1.0 with data from 



15l | and Ass data from pp2pp |16|] . 



which we show in Fig. 19 for the value K = 1.0. For i^ = 0.6 the peaks in each case are 
about 0.6 as high as in Fig. 19. It is clear from this that determination of i^ or a validation 
of the model is not yet possible. (Actually the pp2pp data here is for Ass but there is little 
difference between the two predictions.) A convenient way of parametrizing the strength of 
the double-flip amplitude is through p| 

02(S,O) 



r2{s) 



(6.4) 



2Im0i(s,O)' 

Note that the transverse total cross section asymmetry Acr^ = —2lTar2crtot- ^2 is plotted 
over the fixed target range in Fig. 20 along with la errors. (For K = 0.6 the value of r2 is 
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FIG. 20: Model prediction of r2{s) parameter K = 1.0 with data from [15] . 



at 



about 0.6 times the magnitude.) Its very small value 
far smaller than the errors of the recent experiment 16l] 



s = 2002GeV2 is -0.0005 + 0.0002i, 
. This will be a challenge to observe. 



6-1. The Odderon 



Several years ago Elliot Leader and the author proposed using A^j^ as a means for 



searching for the Odderon pL7|. The idea behind this was that, because the Odderon is 
odd-signature it would be asymptotically real in contrast to the Pomeron 32| . The formula 
for Ann, Eq. 2.4, contains i?e(020i) and there the signal would be enhanced at small t due 
to Coulomb-nuclear interference. This gives a very characteristic signature that we proposed 
using to look for the Odderon. Since at that time we had no information at all about the 
spin flip couplings, our discussion was purely qualitative. Now we have a limited amount of 
information from the preceding discusson and can see if there is any sense to our proposal. 

So we simply add to the expression for 02 a simple Odderon pole at J = 0.96 as a good 
guess based on QCD from Nicolescu [39|, a spin-flip coupling of unknown value Po and zero 
non-flip coupling since there is no sign of it in unpolarized experiments 40|] . We then process 
the new 02 through the absorption model and predict Ann- 

Since we have no idea what the Odderon spin-flip coupling should be, we show in Fig. 21 
the prediction using K = 0.6 for four values of Po, (0,1,2,3), for the three energies where 
there is data. (For K = 1 the prediction is about a factor of 5/3 larger.) From this 
we conclude that we cannot rule out an Odderon with a modest spin-flip coupling, but it 
clearly must be smaller than the normal Regge flip couplings. (Recall that /5£ = 4.67 and 
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Pi = 27.65 for comparison.^ 
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FIG. 21: Prediction for Ann or Assiwith K = 0.6 at (a) pL = 24GeV/c, {b)pL = lOOGeV/c, 
(c) s = 200^GeV with data from 15|,l3a] with Odderon coupling f3o = O(solid),l(short), 2(niediuni) 
and 3 (long) . 



SUMMARY 



This paper serves as a report on the large amount of work done in conjunction with the 
polarized proton program at RHIC. In addition to providing a model which works at the 
10% level for the energy dependence of the the analyzing power of elastic scattering needed 
for proton polarimetry, it also provides some useful information about the spin dependence 
of dominant Regge poles. Most notably, the data indicate that the Pomeron has a small 
but significant spin-flip coupling. These results allow the exploration of the double spin flip 



asymmetry Ann for which some data over a wide energy range are now available [15|, Il6| , 
along with a concrete realization of a proposed Odderon search 17|. Our results are limited 
by the rather large errors on some of the data we have used and we expect that this will 
be remedied in the not too distant future, leading to more precise polarimetry and better 
determination of the Regge spin- flip couplings and so better values for Ann- 
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